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Abstract A model is proposed to study the possible pairing structures of N-boson 
systems with nonzero spin. Analytical solutions have been obtained. The emphasis 
is placed on the spin-structures of ground states with attractive or repulsive pairing 
force, and with or without the action of a magnetic field. A quantity (an analogue of 
the two-body density function) is defined to study the spin-correlation between two 
bosons in N-body systems. The excitation of the system has also been studied. 

1 Introduction 

Historically, the pairing structures played a very important role in fermion systems. 
The most famous example is the Cooper-pair in superconductivity. [I] One more ex- 
ample is the pairing of the valence nucleons in nuclei as described by the "Interacting 
Boson Model" [2,. However, there is no obvious evidence to support the existence of 
paring structures in boson systems. Nonetheless, we can not rule out the possibility. 
For an example, each 52 Cr atom has a spin J = 3. The S-wave scattering lengths a s 
of these atoms depend on the spin-channels s, namely, the total spin of the two inter- 
acting atoms. There are four channels, s = 0, 2, 4, and 6. Accordingly, for low energy 
scattering, the interaction between the i and j atoms can be written as 

Vij = S(r l - r^^gsWij (1) 

where the strength g s = 47rfi 2 a s /M, M is the mass of atom. ^p?- is the projection 
operator of the s— channel. Let x(*) be the spin-state of i, (x(i)x(j))sm = (ij)sm 
be the normalized total spin-state of i and j , where the spins are coupled to s 
and m (the Z-component of s) via the Clebsch-Gordan coefficients. Then, ?Pfj = 
12m \{ij)sm){{ij)sm\- For 52 Cr, it has been determined experimentally that = 
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~7aB, C14 = 58as, and ag = 112as, while ao is undetermined. [3] Since <?4 and <?6 
are strongly repulsive, the atoms in low-lying states would do their best to avoid form- 
ing s — 4 and 6 pairs. Although §2 is attractive, it is weak. Therefore, if go is also 
attractive and strong, the low-lying states might be dominated by the force of the 
s = channel (the pairing force) to form s — pairs. Therefore, in view of having 
so many different kinds of boson systems (the constituent may be an atom, a pair 
of fermions, an exciton, etc.), the assumption that some boson-systems in some cases 
might be dominated by pairing force is assumable. 

In what follows, we shall propose a model to study the effect of the pairing force 
on boson systems. In this model, each particle has a nonzero integral spin J, and they 
are interacting via pairing interaction. In general, if one can obtain analytical solutions 
of a Hamiltonian, then these solutions are valuable because the underlying physics can 
be thereby better understood. Besides, as an exact solution, it can be used to check 
the validity of various approximate methods. It turns out that the above model can be 
solved analytically as shown below. 

Since the experimental realization of the condensation of atoms with nonzero 
spins, [4] the study of the spinor condensates has become a hot topic due to academic 
interest and their potential in application. Traditionally, the mean-field theory is used 
to explain related phenomena. Up to now this theory is successful. Nonetheless, if the 
details of spin-spin correlation are taken into account, this theory is not sufficient due 
to the basic approximation inherent in the theory. On the other hand, the theoretical 
method used in this paper is beyond the mean-field theory. In what follows, approx- 
imation is made only in the spatial space, while the spin-degrees of freedom of the 
iV— body system have been exactly treated. Thus, the following approach would be a 
useful complement to the mean-field theory. The spin-spin correlation might by thereby 
better understood. 



2 Hamiltonian, spectrum, and the eigenstates 

We first consider a N-boson system, each boson has spin J = 3. The cases with J ^= 3 
will be discussed later. The bosons are confined by a potential. The interaction acting 
among the bosons contains a spin-independent term and a spin-dependent term. The 
latter is simply the pairing force, namely, 

Vij =S(ii- Tj)g y$j . (2) 

We are only interested in the case of the low-temperature limit. At this limit all 
particles are condensed into the same spatial state /(r).[S] In other words, the spatial 
state is fixed. Therefore, the spatial degrees of freedom can be integrated. After the 
integration, the Hamiltonian of the model reads 

Hraod = ^2 GoVij (3) 

i<j 

where Go = go J dr\f(r)\ 4 . For this Hamiltonian, only the spin-degrees of freedom 
are considered. Obviously, /(r) depends on the potential of confinement and on the 
interaction. However, since Go will be considered as a parameter, the details of the 
potential and the spin-independent interaction is irrelevant. In what follows, go (and 
therefore Go) is first considered to be negative. 
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For the diagonalization of H mo( i, we introduce the set of symmetrized and normal- 
ized Fock-states jfc) = \n£ , N$, JVf, N$ , N^, N^ 2 , N* 3 ) as basis functions, where 
is the number of atoms with spin component /i, and k denotes the set {NjZ} which 
together determine the basis function. Obviously, ^ — N, the total number of 
particles. Since only the spin-degrees of freedom are involved in H mo( [, the set of Fock- 
states is complete. The related matrix elements read 



(k'\H mod \k) = ^ £(-l)" +, yJV*JV*'[ 5 vfi 5^5 k , k {N% - 1) 

— I — T, ; N k ' — 1 N k ' —1 

+5 vfi 5^NV v {N$-l)5 N l_ 2 - 

+w,oi/^rf:ift i ( 4 ) 

— N k — l,iV* -1 

where = (1 — S^.u). For the label <5 £ , the superscript implies a revised 

set of the set {N$ • • • N* 3 } by reducing both and n!L v by 1, the subscript implies 
a revised set of {N 3 • • • N_ 3 } by reducing iVg by 2. When the two revised sets are 
one-to-one identical, the label is 1, otherwise it is zero. And so on. When N is given, 
the dimension of the matrix is finite. Due to the special structure of the matrix, it can 
be analytically diagonalized. The results are the follows. 

For convenience, N is assumed to be even unless specified. The ground state W g 
has an analytical form as 

SPg = ^^C"j?|fc), (not yet normalized) (5) 

k 

Cl = ^,^ 3 ^,^ 2 V, W * 1 (-l) W3i+Arf ^/((</2)! 2</ 2 ) (6) 

where Nq should be even. On the other hand, we define 

!Pj° ir = P{(12)o(34) • • • (N - 1, AOo} (7) 

where P implies a summation over all the N\ permutations of particles, and the im- 
plication of {ij)o has been given in the previous section. This state can be expanded 
in terms of the Fock-states as 

^ = { -j^(N/2)\2 N l 2 ^N\Y J C 9 k \k) . (8) 

Hence, \pP mr is just the ground state that we have obtained from the diagonalization. 
Thus the structure of the ground state is clear, it is just a product of the s = pairs, 
i.e., a pairing structure. 

We use I Go I as the unit of energy in the follows, then the ground state energy turns 
out to be Eg = — = £n . 

If N is odd, the ground state would be just the above pairing structure together 
with an additional particle, and E g — ej^ + 3/7. 
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It is reminded that the interaction keeps the total spin S and its Z-component M 
to be conserved. Obviously, due to the pairing structure, the ground state has S = 0. 
The first excited level is found to contain the eigen-states having S = 2, 4, and 6. They 
are exactly degenerate. Their analytical forms read 

&y,SM — C%' SM \k), (not yet normalized) (9) 

where the script y implies that the state is a yrast state, and the Clebsch-Gordan 
coefficients have been introduced, and 

F kv = <5L 3 ,L- 3 5L 2 ,L- 2 '5L 1 ,L_ 1 (-l) L3+il N/^!/((io/2)!2 i0 / 2 ) 

y/(L„+ 2)(Lu + l) + ^+l)(L M -„ + l)] (11) 

where L M and TV^; are related as = TV^ — 8^. v — &u,M—v ( A* i s fr° m -3 to 3) and Lq 
should be even. 

Similarly, we define an excited pairing structure by changing a s — pair to a s / 
pair. We found 

P{(12) SM (34) ---(iV-l,iV)o} = C on S x^C^ SM |fc) . (12) 

k 

Thus, it is clear that the Hty SM states have just the excited pairing structures. They 
all have the same eigenenergy E g + 1 disregarding S and N, and the gap is just 1 (in 

|G |, m )■ 

In general, the spectrum of this N— body system contains a finite number of levels, 
each is highly degenerate (except the lowest one) . We introduce a quantum number I 
to denote the levels, therefore I describes the degree of excitation. Each eigenenergy is 
a sum of two terms, it reads 

E I = e N + e I . (13) 

Where £jy = ~ jv ^ +J -* depends roughly on the number of pairs (~ TV 2 ) and is the 
same for all the eigen-states, while El = I -^T L depends on /, J = 0, 2, 4, • • • N (if 
N is even, in this case the first excited level has 7 = 2 and therefore ej = 1), or 
I = 1,3, •••N (if N is odd). So, when N is even (odd), totally there are N/2 + 1 
((iV+l)/2) energy levels. When / is larger, the level will be not only higher, but also 
have a larger degeneracy as shown later. When we compare the spectra of two systems 
distinct in TV, their ej are one-to-one identical (except the range of 7). Thus, although 
the spectrum of the one with more particles would contain more levels, the lower parts 
of their spectra are exactly the same except the shift caused by e^. This is a noticeable 
point. The spacing of adjacent levels is just — Ej — (21 + 7)/7. Thus, when I 

goes up, the spacing becomes larger and larger. Note that, since the set of Fock-states 
is complete, the above spectrum is exact for the model Hamiltonian. 

For all the states of the I-level, I atoms would be excited from the s = pairs. 
Therefore, I has a similar implication as the seniority which has been introduced for 
describing the excitation of pairing structures of fermion systems long ago. 7,8 Fur- 
thermore, the maximum S of the states belonging to an I-level is obvious 3/ contributed 
by the I excited atoms. 
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In general, the eigen-states can be classified according to S, M, and 7 . It is possible 
that there are more than one states having the same set (ISM) (e.g., for the 1 = 4 
level, there are two independent S = 4 states). In this case an extra quantum number 
is necessary. Therefore, an eigenstate can be labeled as ^i^s,M,iy where i denotes other 
good quantum numbers, if necessary. For examples, the above three yrast states ^ Vt SM 
have 1 = 2, therefore they can also be denoted as l?2,SM- 

The degeneracy of an I- level is stated as follows. Let {j^} be a set of seven non- 
negative integers (/i is from -3 to 3). There would be a finite number of sets satisfy 
both = / and aj^ = 0. The number of these sets is denoted by Oj. Then, 

one can prove that the degeneracy of the I-th level di = Qj — (do = d\ = 1). 

This number increases very fast with I, e.g., do = 1, <fe = 3, and d 4 = 14. Note that 
the degeneracy caused by M has not yet been taken into account. 

For higher excitation, a stste from the 1 = 4 level is given as an example as 

^400 = cons x (3P{[(12) 2 (34) 2 ] (56)o ■ ■ • } - ^P{[(12) 4 (34) 4 ]o(56) ■ ■ ■ }) (14) 

where four particles are excited from the s = pairs. The structures of higher states 
will become more and more complicated because more particles are excited. We are 
not going to the details of them. 

When Go is positive, the above analytical solutions hold exactly. However, the 
spectrum would be reversed. The ground states would be highly degenerate, all of 
them belong to the I = N level, wherein all s = pairs are rejected. In particular, 
there is a great energy gap E gap = £?/v-2 — = (2N + 3)/7. A larger N leads to a 
larger gap. Thus, when N is large, the spin-degrees of freedom of the group of ground 
states are difficult to be excited. This is also a noticeable point. 

When J / 3, one can also use the Fock-states as basis functions. After the di- 
agonalization, the eigenenergies are also a sum of two terms, namely, Ej = en + ej. 
However, = — ~^pqrn ~ > wnne £ I = ~^T7~fT7^' e ig ens tates have the same 
structures as before, and therefore can also be classified as ^iSMi- in particular, the 
1 = level has all particles in s = pairs, while the 1 = 2 level has two particles 
excited from the s = pairs, and so on. 

3 Ground states under a magnetic field 

We go back again to the cases with J = 3 and an even N and a negative Go- When 
a magnetic field B lying along the Z-axis is applied, the model Hamiltonian reads (in 

IGol) 

H' mod = -Y,y>%-lS z (15) 

i<j 

where Sz is the operator of the total spin along the Z-axis, 7 = g^igB /\Gq\, fis is the 
Bohr magneton (incidentally, if the atom is 52 Cr with J = 3, g = 2). The eigen-states 
of H' mod are exactly the same as H mod , however the eigen-energies would contain an 
additional term due to the Stern-Gerlach splitting as 

E'ism = £N + £1 - lM . (16) 

We are interested in how the structure of the ground state would vary with B. For 
the group of states belonging to the 7— level, the largest M possessed by these states 
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is 37. Thus the state with S — M = 3/ of each J— level is a candidate of the ground 
state with the energy E\ 3I 3I . Among the candidates, when B (or 7) is fixed, if I = I 
leads to the minimum of the group Ej 3/ 3/ , then the ground state would have I — I - 
Obviously, I depends on B. It was found that, when B is small so that 7 < 7/42, 
the ground state has I = 0. Starting from 7/42, let the range of 7 be divided into 
segments, each has a length 2/21. Then, when 7 increases, I will increase step by step 
from a segment to the next segment. Each step I will jump by 2. The general relation 
of To and 7 is: 

(2J + 3)/42<7< (2 1 + 7) /42 . (17) 

Since I Q can not exceed N, the last segment is (2N + 3)/42 < 7 < 00. In this 
extensive segment, I = N, namely, the increase of the magnetic field will eventually 
excite all particles from the s — pairs. 

In each segment the associated ground state energy 75 ffr (7) = E' Io 3J(> 3I ^, and 
-j^Egri'y) = — 3/ within the segment. When 7 increases, since I will increase step by 
step, the slope becomes more and more negative. 

When 7 < (2N + 3)/42, we have approximately 7 w (27 c + 5)/42. Thus we have 
Egr^) « — £at — (427 — 5) 2 /56. It implies that, when 7 increases, the decline of E gr (^) 
against 7 is nearly parabolic. Of course, when 7 > (2N + 3)/42, the decline is only 
contributed by the last term of eq.(16), i.e., —3~/N, and therefore is linear. 

The analytical form of the wave function of the ground state can also be obtained. 
It reads in terms of the Fock-states as (not yet normalized) 

-VM = *i e ,3i e ,3i = E*iVs*-Io,iv^*iv^iv^^^iv^(-l)^ +w V^o w 3 ! /Jv*3 ! 
k 

/((7V fe /2)! 2< /2 ) \k) 
= cons x P{[/T/' 1 X3(t)](/o + 1, Io + 2)o • ■ ■ (N - 1, iV) } . (18) 

Thus, when a magnetic field is applied, the ground state would have a portion of 
particles excited and become fully polarized, while the others remain in the s = pairs. 
Thus the implication of I is clear, it is just the number of polarized particles. 

From &gr("i) it is straight forward to calculate the probability of a particle in 
/i— component P M = (</>| 7V M /AT, where Np is the operator of the number of particles 
in u, and — 1- Obviously, if P3 = 1 (or P_3 = 1) while the other P M = , 

the system is fully polarized. Therefore, from one can understand the situation of 
magnetization of a state. For fl*gr(7), from (18) we have 

D I (Io + N + 6)+N 
^3 = 



P-3 = 



N(7 + 2I ) 

(Iq + 1)(N-Iq) 

N(7 + 2I ) 



P ^wrkry (if ^ ±3) (19) 

Examples of against 7 are given in Fig.l. When 7 = 0, all P^ — 1/7. In this 
case all particles are in s = pairs. Accordingly, the total spin-state has 5 = and 
therefore is isotropic. When 7 increases, the variation of P^ can be roughly divided 
into three stages. In the first stage both P3 and P_3 increases very fast while the other 
Pfi decreases very fast as shown in the figure. When 7 = 93-/V/(21(2iV + 100)) = 71, 
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Fig. 1 P M (the probability of an atom in fi component) of the ground states of H' , against 
7 with N = 100 (a) and 10000 (b,c). P3 is in solid line, P_3 in dash line, and all other are 
in dash-dot-dot line (they overlap exactly), (c) is just the left end of (b). 



Pfi (/•* 7^ ±3) are found to be equal to 0.01. We use the range to 71 to define the first 
stage. The second stage is from 71 to 72 = (2N + 3)/42 (say, when N = 10 4 , 71 =2.2 
and 72 = 476.3 as shown in Fig.lc and lb). In this stage P3 increases (PL 3 decreases) 
nearly linearly until equal to 1 (0), while the other P^ remain w 0. The third stage 
(7 > 72) is characterized by P3 = 1, i.e., the system is fully polarized. 

Incidentally, the structure of the ground states of the 52 Cr condensates under a 
magnetic field has already been studied by a number of authors using the mean field 
theory. [9]?,?] Obviously, the spinors from the mean field theory can be compared with 
Pfi. When go is very negative, two phases "polar" and "ferromagnetic" have been found. 
The associated spinors are (cos 9, 0, 0, 0, 0, 0, sin 9) and (1,0,0,0,0,0,0), respectively. Ob- 
viously, the former (latter) corresponds to the above second (third) stage of P^. From 
the correspondence we know that the underlying physics of the "polar phase" is just 



8 



a kind of structure wherein a portion of particles are fully polarized while all other 
particles are in s = pairs. 

Using the concept of the fractional parentage coefficients one can extract any pair 
of particles, say i and j, from a spin-state 9 — ~}2 k C\\k} as 



* = E \ / < ( 1 > \ - 1 ./ ) [ v< E c * \/ x j } i • • • ■ N Z - 2 > • • • ) 



+^E^J^)|- ,^-1,...)] (20) 

where the two Fock-states written explicitly at the right are for (JV — 2)-body systems. 
The former originates from \k) by changing N% to iV^ — 2, while the latter by changing 
both and AT* to iV* — 1 and N„ — 1, respectively. From this expansion it is straight 
forward to calculate how the magnetization m of a selected pair would be. Let /i and ^ 
denote the spin-components of the two selected atoms i and j, and m = // + ia Then, 
the probability of the pair in m reads 

9™-^W,">IVL c * jv(jv-i) m "y fe iV(iV-i) J ' ( } 

In addition to Q m is a measure of spin-correlation (corresponding to the two- 
body correlation function in spatial space), and has no counterpart in mean field theory. 
Qm of the ground states of H' mod against 7 are plotted in Fig. 2. Among all the thirteen 
Qm (m is from -6 to 6), Qo is much larger than the others if 7 < 7/42 (in this case 
all particles are in s — pairs). Thus, for every pair of particles, the probability of 
having their spins anti-parallel is the largest. The next largest is the group Q\ — Q2 = 
Q-i — Q-2, while the smallest is the group Qe ~ Q5 = Qs = Q-5- However, when 
7 increases, Qq will increase nearly linearly. Qq and Q-q increase firstly with 7 but 
eventually decrease with 7 as plotted. Other Q m are very small. When 7 > 72, Qe = 1 
implying a full polarization as expected. 



4 Summary 

A model is proposed to study the possible pairing structures of N-boson systems with 
nonzero spin. Analytical solutions have been obtained. The following points are men- 
tioned 

When the pairing force is negative, the ground state has all particles forming the 
s — pairs. The excited states are described by good quantum numbers S, M, and I. 
The latter is just the number of particles not in the s = pairs. The state has a larger 
/ is higher, and the states have the same I are degenerate. When a magnetic field B 
is applied, the ground state will have a portion of particles fully polarized, while the 
other particles remain in the s = pairs. 

When the pairing force is positive, the ground state is highly degenerate with a gap 
proportional to N. Therefore, when N is large, the spin-degrees of freedom are difficult 
to be excited. 

Two quantities and Q m have been defined and studied. They are the spin-space 
analogues of the one- and two-body densities of coordinate space. In particular, Q m 
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can help us to understand better the spin-correlation and has no counterpart in mean 
field theory. Therefore, these quantities are in general useful for understanding the spin- 
structures of various systems. For an example, the underlying structure of the "polar 
phase" found in |9j?,?] is thereby clarified as a mixture of fully polarized particles and 
s — pairs. 

In general, the results from this model system might be helpful for the understand- 
ing of some realistic systems. 
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